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Abstract. We will discuss recent results on the thermodynamics of QCD in the presence of light 
dynamical quark degrees of freedom. In particular, we will concentrate on an analysis of the flavour 
, and quark mass dependence of the QCD phase diagram, the equation of state and the transition 

temperature. Moreover, we present recent results on the heavy quark free energy. 

<n : 

INTRODUCTION 

o : 

q ! Lattice calculations have provided a rather detailed picture of the thermodynamics of 

Q\ ; gluonic matter. We know that a phase transition exists, which separates a confining low 

temperature phase (glueball gas) from a deconfined gluon gas phase at high temperature. 
The transition between these two phases is first order, as originally predicted by Svetit- 
"^3 | sky and Yaffe [1]. Although the low temperature phase only consists of rather heavy 

glueballs, the lightest of which has a mass of about 1.8 GeV [2], the phase transition 
^ ■ temperature turns out to be rather small, T c = 0.637(5)^/o ~ 270 MeV. This suggests 

<D | that T c does not sensitively depend on the mass of the lightest excitations in the confining 

^ ■ phase. In fact, it agrees quite well with values extracted from resonance gas models based 

on the excitation spectrum of a gluonic string [3]. Bulk thermodynamic observables like 
the energy density and pressure change rapidly at T c and asymptotically approach the 
ideal gas limit. Their rapid rise signals the liberation of many new light degrees of free- 
dom at T c . However, even at temperatures a few times T c one still observes significant 
deviations from the asymptotic ideal gas behaviour. This suggests that also in the high 
temperature phase non-perturbative effects play an important role, which give rise to 
large screening masses and quasi-particle excitations. 

Gluonic matter is described by an SU (3) gauge theory, which is obtained as the infinite 
quark mass limit of QCD. In this limit the entire fermionic sector of QCD decouples 
and does not contribute to the thermodynamics; quarks only serve as static sources 
(quenched QCD). This allows, for instance, a study of thermal modifications (screening) 
of the forces between external charges, which shows that the linear confining potential 
weakens with increasing temperature; the string tension decreases and vanishes at T c . 

Some of these basic aspects of gluonic thermodynamics clearly will change drasti- 
cally in the presence of light dynamical quark degrees of freedom. The asymptotic high 
temperature limit for bulk thermodynamic observables like the energy density or pres- 
sure rises with increasing number of light degrees of freedom. Absolute confinement is 
lost in the presence of quarks with finite mass already at low temperatures; the heavy 
quark free energy will no longer diverge when one tries to separate a static quark anti- 
quark pair. Moreover, the order of the QCD phase transition and even its very existence 



will crucially depend on the number of light degrees of freedom and their masses. 

During recent years these basic qualitative changes in the thermodynamics of QCD, 
which result from the presence of light quark degrees of freedom, have been observed 
in lattice calculations [4, 5]. We will discuss here the current status of their quantita- 
tive analysis. This will make clear that unlike in the quenched sector of QCD we did 
not yet reach a similarly detailed quantitative understanding of the relevant parameters 
that control the critical behaviour of QCD with a realistic spectrum of quark masses. 
However, thermodynamic calculations on the lattice steadily improve. This partly is due 
to the rapid development of computer technology. However equally important has been 
and still is the development of improved discretization schemes, i.e. improved actions. 
This is of particular relevance for thermodynamic calculations which are not only sensi- 
tive to the long distance physics at the phase transition but also probe properties at short 
distances in calculations of e.g. the energy density or the heavy quark potential [6]. 

In the following we will discuss some recent results on the flavour and quark mass 
dependence of QCD thermodynamics. We will not go into any details on the lattice 
formulation of QCD at finite temperature, which have been discussed elsewhere [7]. 



THE QCD PHASE DIAGRAM 

The quark mass and flavour dependence of the QCD phase transition at finite temper- 
ature and vanishing baryon number density has been studied extensively in lattice cal- 
culations. The basic qualitative and quantitative features expected from universality ar- 
guments [1, 8] on the one hand and phenomenological considerations on the other hand 
have been reproduced by these calculations. The transition is found to be first order in 
the limit of infinitely heavy quarks as well as in the chiral limit of 3-flavour QCD. In 
the case of 2-flavour QCD the transition is found to be continuous. The current status of 
the analysis of universal properties in the chiral limit, however, is not really satisfying 
[9]. The demonstration that for 2-flavour QCD the transition belongs to the universality 
class of 3-d, 0(4) symmetric spin models still is ambiguous. 

The regions of first order transitions are separated from a broad crossover region 
by lines of second order phase transitions. These lines are expected to belong to the 
universality class of the 3-d Ising model [10], which is quite remarkable as the global 
Z(2) symmetry which gets restored at these transitions is not a symmetry of the QCD 
Lagrangian. This makes it obvious that the bare couplings appearing in the QCD La- 
grangian cannot be the relevant scaling fields that control the critical behaviour at these 
lines of second order transitions. It therefore also is not obvious a priori what is the 
correct order parameter for these transitions. The energy-like and magnetization-like 
operators of the effective Hamiltonian that controls the critical behaviour at this critical 
point will be linear combinations of the basic fields appearing in the QCD Lagrangian 
and the relevant temperature-like and ordering-field like couplings will be linear combi- 
nations of the bare parameters of the QCD Lagrangian, i.e. the gauge coupling (3 = 6/g 2 
and the quark masses m q . 

Understanding the critical behaviour and the relevant scaling fields in the vicinity 
of the chiral critical line at small values of the quark masses is, of course, interesting 




in its own rights. However, it eventually will also be of importance for a discussion 
of the physics in the vicinity of the critical endpoint which is expected to exist in the 
temperature-density phase diagram [11]. In an extended phase diagram, which also 
includes the dependence on the chemical potential, these critical points lie in the same 
critical surface of second order transitions. 

The critical point separating the first order from the crossover region has recently been 
studied in some detail for the case of three degenerate quark mass flavours [12]. Through 
an analysis of cumulants of the chiral condensate, (\jA|/), as well as joint probability 
distributions for the chiral condensate and the gluonic part of the QCD action it could be 
shown that the transition indeed belongs to the universality class of the 3-d Ising model 
and that the corresponding order parameter can be constructed as a linear combination 
of the chiral condensate and the gluonic action. In Figure 1 we show the result of a 
calculation of the fourth cumulant of the chiral condensate (Binder cumulant), 
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The cumulant has been calculated for different values of the quark mass at the pseudo- 
critical couplings. Up to finite volume corrections the cumulants obtained on different 
size lattices will cross at the critical quark mass corresponding to the second order phase 
transition point. The value of B4 at this point is universal and unambiguously identifies 
the transition as an Ising-like transition. 

In order to judge whether QCD with a realistic spectrum of light u,d quarks and a 
heavier strange quark lies in the first order or crossover region of the phase diagram it 
is necessary to determine the location of the critical line quantitatively. The calculations 
performed so far with 3 degenerate quark masses and different fermion actions [12] show 
that the critical mass parameter, e.g. the value of the pseudo-scalar meson mass at the 



chiral critical point, is rather sensitive to cut-off effects. While the calculation performed 
with the standard staggered fermion action led to a critical value m ps ~ 300 MeV, 
calculations with an improved staggered fermion action gave m ps ~ 200 MeV. These 
calculations thus need further confirmation through studies closer to the continuum 
limit. Nonetheless, the current estimates consistently yield rather small values for the 
pseudo-scalar meson mass at the critical endpoint in 3-flavour QCD, which makes it 
quite unlikely that the physical point in the phase diagram, corresponding to a realistic 
quark mass spectrum with two light u,d, and a heavier strange quark, would lie in the first 
order region. This would also be in agreement with an earlier estimate of the Columbia 
group [13]. 

Our current understanding of the QCD phase diagram of 3-flavour QCD at vanishing 
baryon number density is summarized in Figure 2. 
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FIGURE 2. The QCD phase diagram of 3-flavour QCD with degenerate (u,d)-quark masses and a 
strange quark mass m s . 



THE TRANSITION TEMPERATURE 

It should be clear from the previous discussion of the phase diagram that for a large range 
of quark masses the transition to the high temperature plasma phase is not a phase transi- 
tion, i.e. the transition does not correspond to singularities in the QCD partition function. 
Nonetheless, also for these quark masses the transition occurs in a narrow temperature 
interval and thus is well localized. This transition region is characterized by peaks in 
susceptibilities, e.g. the chiral susceptibility % m or the Polyakov-loop susceptibility %l, 

5U = ^<W> , XL = N 3 a [(L 2 )-(L) 2 ) . (2) 



Here 




£TrL f 



(3) 



denotes the spatial average over Polyakov-loops 1 , Lf, defined at the spatial sites x 
of a lattice of size iV| x A^. A calculation of the pseudo-critical temperature, T c = 
l/N x fl(Pp C ) still requires the determination of the lattice spacing a((3 pc ) at the pseudo- 
critical couplings $ pc which in turn are defined through the location of the susceptibility 
peaks. The lattice spacing can be determined through the calculation of an independent 
physical observable. In order to quote T c in physical units, i.e. MeV, we need a physical 
observable which does not crucially depend on the values of the quark masses. For 
physical values of the quark masses it would, of course, be most appropriate to use 
a hadron mass, e.g. the rho-meson mass, to set the physical scale for T c . The hadron 
masses, however, are themselves strongly dependent on the quark mass values; their 
masses diverge in the infinite quark mass limit. Nonetheless, we can assign a physical 
value to the transition temperature in this limit. In quenched QCD the natural scale for T c 
is the square root of the string tension, In fact, the string tension as well as quenched 
hadron masses 2 seem to describe the experimentally known QCD spectrum reasonably 
well already in the infinite quark mass limit. These quantities thus seem to have a weak 
quark mass dependence and are suitable to set the scale for T c at arbitrary values of m q . 
This situation is illustrated in Figure 3. 

On the left hand side of Figure 3 we show the transition temperature in units of the 
vector meson mass my versus the ratio m ps /my, which in the chiral limit is proportional 
to the square root of the quark mass. Here the drop in T c /my observed with increasing 
quark mass mainly is due to the quark mass dependence of my and does not reflect the 
quark mass dependence of T c . This figure shows, however, that calculations based on 
different discretization schemes for the gauge and fermion actions do yield consistent 
results for T c and its dependence on m q . In the chiral limit one finds for the critical 
temperature in 2 and 3-flavour QCD 



3 — flavour QCD : T c = (154 ± 8) MeV, improved staggered fermions [15] 

Here m p has been used to set the scale for T c . We note that all the results presented in 
Figure 3 have been obtained on lattices with a rather small temporal extent (N z = 4). The 
lattice spacing is therefore still quite large in the vicinity of T c , i.e. a ~ 0.3 fm. Moreover, 



The Polyakov-loop is a purely time-like Wilson loop which is closed due to the periodicity of the 
lattice in temporal direction. Its definition as well as further details on the lattice formulation of QCD 
thermodynamics may, for instance, be found in [7] . 

2 In the calculation of quenched hadron masses the contribution of dynamical light sea quarks is sup- 
pressed. Only the valence quark contributions and interactions with the gluonic vacuum are taken into 
account. 



2 - flavour QCD : 




(171 ±4) MeV, clover-improved Wilson fermions [14] 
( 173 ± 8) MeV, improved staggered fermions [15] 




FIGURE 3. Transition temperatures in units of my (left) and in units of \[g (right). The left hand 
figure shows results for 2-fiavour QCD obtained with different gauge and fermion actions. The right hand 
figure shows the transition temperature in 2 (filled squares) and 3 (circles) flavour QCD obtained with 
improved staggered fermions (p4-action). Also shown are results for 2-flavour QCD obtained with the 
standard staggered fermion action (triangles). The dashed band indicates the uncertainty on T c / ^/a in the 
quenched limit. The straight line is the fit given in Eq. 4. The vertical lines correspond to physical values 
of the hadron masses and a string tension of 425 MeV. 



there are uncertainties involved in the ansatz used to extrapolate to the chiral limit. We 
estimate that the systematic error on the value of T c /m p still is of similar magnitude as 
the purely statistical error quoted above. The critical temperature thus is at present only 
known with an error of about 10%. 

On the right hand side of Figure 3 we give T c in units of ^/o. This reflects the expected 
behaviour; with increasing values of the quark mass the hadrons become heavier and it 
becomes more difficult to create a dense hadronic system that could undergo a transition 
to a quark-gluon plasma phase; the transition temperature thus increases. It is, however, 
striking that the quark mass dependence of T c is so weak. The straight line shown in 
Figure 3 is a fit to the 3-flavour data, which gave 

(T c /Vo) x = (T c /y/a) + 0.04(1)* with x = m ps /Vo . (4) 

This suggests that also heavier resonances, which are little affected by changes of the 
quark mass, play an important role for building up the critical density needed to trigger 
the QCD phase transition. Moreover, we note that for pseudo-scalar masses larger than 
about 6\f<5 ~ 2.5 GeV the transition temperature agrees with the value found in the 
pure gauge theory, T c ~ 0.637^/a ~ 270 MeV. In this heavy quark mass regime all the 
hadronic states are heavier than typical glueball masses and thus decouple from the 
thermodynamics, which only is controlled by gluonic degrees of freedom. In this heavy 
mass regime the transition also changes again from a crossover to a first order phase 
transition. 



THE QCD EQUATION OF STATE 



The temperature dependence of bulk thermodynamic observables like the energy den- 
sity (£) and pressure (p) has been analyzed in detail in the pure gauge sector. These 
calculations show that e/T 4 as well as p/T 4 rapidly increase above T c . However, even at 
T —4T C the Stefan-Boltzmann limit is not yet reached [16]. Deviations stay at the level 
of 15%, which is too much to be understood in terms of perturbative corrections [17]. 
Even at these high temperatures non-perturbative effects seem to play an important role 
for the thermodynamics. This finds further support from studies of the heavy quark free 
energy [18], the spatial string tension [19] as well as the gluon propagator in fixed gauges 
[20]. The analysis of the temperature dependence of these observables also suggests that 
the temperature dependent running coupling constant remains large in the plasma phase 
and electric and magnetic screening masses differ significantly from perturbative results 
even at temperatures which are several orders of magnitude larger than T c [21]. In view 
of this it may even be surprising that models with weakly interacting quasi-particles [22] 
as well as calculations based on self-consistent resummation schemes [23] do provide a 
quite satisfactory description of bulk thermodynamics down to temperatures a few times 



At least at high temperature, where the energy density and pressure are expected 
to approach the free gas limit, the QCD equation of state will strongly depend on the 
number of light partonic degrees of freedom. Already for two massless quark flavours 
the Stefan-Boltzmann constant increases by more than a factor two relative to the pure 
gauge theory. For ^/--flavour QCD one has, 



As the influence of non-zero quark masses on bulk thermodynamic observables will 
be exponentially suppressed at high temperature, Eq. 5 also gives the dominant high 
temperature behaviour for massive quark even when the masses are of the order of the 
(critical) temperature. In fact, this has been observed in lattice calculations with fairly 
large quark masses [24]. 

The overall pattern of the temperature dependence of bulk thermodynamic observ- 
ables in QCD with 2 and 3 light quark flavours is very similar to the case of the pure 
gauge theory. With increasing number of light degrees of freedom the critical tempera- 
ture shifts to smaller values and the asymptotic high temperature limit for p/T 4 becomes 
larger. This is seen in the left hand part of Figure 4. However, after rescaling the ther- 
modynamic observables with the corresponding Stefan-Boltzmann constants they look 
quite alike in units of T /T c . Of course, they still differ in details. In particular, it is ap- 
parent from the right hand part of Figure 4 that at T c the rescaled pressure of QCD with 
light quarks is significantly larger than in the pure gauge theory. This also is the case for 
the energy density which is found to be £ C /T 4 ~ 6 in QCD with light quarks [24, 25], 
while it is only £ C /T 4 ~ 1 in the SU(3) gauge theory 3 . Much of this factor 6 difference in 



J The energy density is discontinuous in the case of the pure gauge theory. The latent heat is found to be 
Ae/r e 4 = 1 .40(9) [26] with e/7; 4 ~ 2 in the high temperature phase. 
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FIGURE 4. The pressure in QCD with «/ = 0, 2 and 3 light quarks as well as two light and a heavier 
(strange) quark. For rif ^0 calculations have been performed on a iV t = 4 lattice using improved gauge 
and staggered fermion actions [15]. In the case of the SU(3) pure gauge theory the continuum extrapolated 
result is shown [16]. Arrows indicate the ideal gas pressure p SB as given in Eq. 5. 



e c /r c 4 , however, seems to arise from the difference in T c between QCD with light quarks 
and the purely gluonic theory. The critical energy densities turn out to be quite similar. 
Unfortunately, the current error on T c , which is about 10%, amplifies in the calculation 
of the energy density, which makes e c still badly determined in lattice calculations, 

g c = (0.3-1.3) GeV/fm 3 . (6) 



THE HEAVY QUARK FREE ENERGY 

The confining properties of the thermal heat bath generated by quarks and gluons can 
be probed by analyzing the response of the medium to the insertion of static sources. 
Static quark sources are described by the Polyakov-loop L%. In the pure gauge limit 
the expectation value of its spatial average, (L) (Eq. 3), is an order parameter for the 
deconfinement transition, 

= , T < T c 



< l h> ; t>£ ■ ™ 

This reflects the long distance behaviour of the correlation function for static quark anti- 
quark sources, 

Q -F{r,T)/T = ( TrL() TrL t^ ^ r = pe| ? ( 8 ) 

which approaches the cluster value | (L) | 2 in the limit of infinite separation between the 
quark anti-quark pair. In the confined phase this correlation function vanishes for r — > °° 
which signals that the free energy needed to separate the two sources is infinite. This is 
no longer the case if dynamical quarks with a finite mass are contributing to the thermal 
heat bath. The static sources can then be screened by quarks and anti-quarks present in 
the heat bath and even in the zero temperature limit this becomes possible through the 
generation of ^^-pairs from the vacuum. The free energy needed to separate the static 
sources thus will stay finite at all temperatures. 




FIGURE 5. The heavy quark free energy in units of the square root of the string tension versus the 
quark anti-quark separation (left) and the estimate for the dissociation energy defined in Eq. 9 (right). The 
band shows the normalized Cornell potential V(r) = — a/r + a r with a = 0.25 ±0.05. The lattice data 
have been normalized to this potential at the shortest distance available, i.e. at rT = 1 /4. 



In Figure 5 we show the free energy of a static quark anti-quark pair calculated 
in 3-flavour QCD at various temperatures and for different separation of the static 
sources. At low temperatures, T<0J T c , the heavy quark free energy coincides with 
the confining Cornell-type potential, V(r) = — a/r + o r with a = 0.25 ±0.05, up to 
distances r ~ 1.5/ \/o ~ 0.7 fm. With increasing temperature F(r,T), however, gets 
screened earlier and at T c it starts to deviate from the zero temperature potential already 
at distances r ~ 0.3 fm. Moreover, it becomes much easier to separate the heavy quark 
sources. As an estimate for the dissociation energy we show in the right hand part of 
Figure 5 the difference in free energy of a qq-pak at infinite separation and a qq-pak at 
distance r q = y/a/a, 

AF=\imF(r,T)-F(Ja/o) . (9) 

r — >°o 

The rapid decrease of AF close to T c clearly will have consequences for the formation 
and existence of heavy quark bound states not only in the high temperature phase but 
also in the vicinity of T c . Already at T ~ 0.9 T c the free energy difference for a heavy 
quark pair separated by a distance similar to the 7/\|/ radius (r v ~ 0.2 fm) and a ^-pair 
at infinite separation is only 500 MeV, which is compatible with the average thermal 
energy of a gluon (~ 3 T c ). The cc-bound states are thus expected to dissolve already 
close to T c , maybe even below T c [27]. 
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